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We develop a Boltzmann-Langevin description of Coulomb drag effect in clean double-layer sys¬ 
tems with large interlayer separation d as compared to the average interelectron distance Af- 
Coulomb drag arises from density fluctuations with spatial scales of order d. At low tempera¬ 
tures, their characteristic frequencies exceed the intralayer equilibration rate of the electron liquid, 
and Coulomb drag may be treated in the collisionless approximation. As temperature is raised, the 
electron mean free path becomes short due to electron-electron scattering. This leads to local equi¬ 
libration of electron liquid, and consequently drag is determined by hydrodynamic density modes. 

Our theory applies to both the collisionless and the hydrodynamic regimes, and it enables us to 
describe the crossover between them. We find that drag resistivity exhibits a nonmonotonic tem¬ 
perature dependence with multiple crossovers at distinct energy scales. At the lowest temperatures. 
Coulomb drag is dominated by the particle-hole continuum, whereas at higher temperatures of the 
collision-dominated regime it is governed by the plasmon modes. We observe that fast intralayer 
equilibration mediated by electron-electron collisions ultimately renders a stronger drag effect. 

PACS numbers: 71.27.+a, 72.10.-d, 73.40.Ei, 73.63.Hs 


I. INTRODUCTION 


Since the pioneering papers by PogrebinskiPl and 
Pricetheir idea of Coulomb drag has evolved into the 
independent and fruitful field of research in condensed- 
matter physics. The initial motivation was to propose an 
experiment that would provide a direct measure of the 
rate of electron-electron collisions not masked by other 
competing relaxation channels, such as electron-impurity 
or electron-phonon collisions. This became possible and 
was realized in the electrically isolated but interactively 
coupled double-layer systemsWhen current is driven 
through one layer while the voltage across the other layer 
is measured, the resulting nonlocal drag resistivity is a di¬ 
rect probe of the rate of momentum transfer between the 
two layers via electron-electron scattering. In many prac¬ 
tical cases, Coulomb drag measurements provide incred¬ 
ibly sensitive tools for revealing electronic correlations, 
which are not readily accessible via more standard trans¬ 
port experiments in single-layer samples. Recent results 
on drag resistivity in double-layer heterostructures made 
of graphen^^KSl triggered a flood of theoretical worlPI®^ 
including new proposals for the mechanism of this phe¬ 
nomenon in the case of tightly nested layers 

Coulomb drag resistivity is extremely sensitive to 
temperature T, magnetic field B, interlayer spacing d, in¬ 
tralayer density n or density mismatch between the lay¬ 
ers, and intralayer mean free path i, which can be domi¬ 
nated by either impurity scattering in the disordered case 
or by electron-electron collisions in clean systems. The 
strength of correlation effects can be conveniently de¬ 
scribed by the electron gas parameter = (7rna^)“^/^, 
where as = ejme^ is the effective Bohr radius in the host 
material of a quantum well with dielectric constant e. 


In initial measurement both the characteristic de¬ 
pendence of the drag resistivity on the various parameters 
and the magnitude of the effect w ere in very good agree¬ 
ment with theoretical predict ions In particular, 

at relatively low temperatures and in the clean limit, 
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where kp and Ep are Fermi momentum and energy, re¬ 
spectively, K Vskp IS the inverse Thomas-Fermi screen¬ 
ing radius, and pg = 27r/e^ is quantum of resistance 
(throughout the paper, h = ks = !)• These early exper¬ 
iments were performed on quantum wells with the typi¬ 
cal electron density n ^ lO^^crn”^ and layer separation 
d ^ 250A. That translates into Ep ^ 60K, kpd^ 1 and 
corresponds to the weakly interacting limit ^ 1, which 
was explored in theoretical calculations. 

Subsequent detailed experimental investigations of 
drag resistivity po sed t heoretical challenges. For exam¬ 
ple, measurementon samples with densities and 
mobilities comparable to that of early measurements, 
but with much higher interlayer separations (up to d ^ 
5000A), revealed that polT^ ceases to be a constant but 
acquires a peculiar temperature dependence, while the 
overall magnitude of the drag resistivity significantly ex¬ 
ceeds the expected value. Such a dramatic disagreement 
between experimental findings and the theoretical predic¬ 
tions of Eq. 0 was attributed to an additional contri¬ 
bution to drag effect due to virtual-phonon exchange.l^ 
It was also observed that at higher temperatures, drag 
resistivity becomes a nonmonotonic function of T. It ex¬ 
hibits a well-pronounced peak at T ^ E p followed by 
a rapid decay at higher temperaturesThe difficulty 
of explaining this feature within the phonon drag mech- 
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anism prompted cons iderat ion of the plasmon contribu¬ 
tion to Coulomb dragP^^^ which yields an enhancement 
of drag at temperatures of the order of the characteris¬ 
tic plasmon energies. Further challenges to the theory of 
Coulomb drag were posed by experiments in samples with 
low carrier density, n ^ lQ^c m~^, w here inter-electron in¬ 
teractions are strong, ^ il 39 | 40 | such double-layers, 
even at low temperatures, T ^ the drag resistance is 
one to two orders of magnitude larger than expected on 
the basis of a simple extrapolation of the small results. 
Furthermore, the power exponent in the temperature de¬ 
pendence deviates from being simply quadratic, drag has 
unexpected behavior on the density mismatch, and the 
system has an anomalous response to a magnetic field. 

For electron liquids with ^ 1 there exists a wide 
temperature interval, Ep < T < VgEp, in which the 
liquid is strongly correlated but is not quantum degen¬ 
erate, and hence it may not be described by the Landau 
Fermi-liquid theory. Within this interval, one can fur¬ 
ther distinguish b etwee n the classical, ujp < T < VgEp, 
and semiquantum JMliH regimes (here 

ujp rsu Ep-s/r^ is the Debye frequency in the liquid). The¬ 
ory of Coulomb drag in this temperature interval has not 
been developed. A detailed consideration of Coulomb 
drag based on an extrapolation of Fermi-liquid-based for¬ 
mulas to the region where > 1 has been carried out 
in Ref. [43] in an attempt to address the data of Ref. (391 
however such analysis can be qualitatively justified only 
for the temperature range T < Ep/(kpd). Another ap¬ 
proach to the theory of Coulomb drag in systems with 
Vs ^ I was developed in the hydrodynamic regime,!^ 
where the mean free path due to intralayer electron- 
electron collisions becomes shorter than other relevant 
length scales. In this case, the drag resistivity can be 
expressed in terms of the viscosity and thermal conduc¬ 
tivity of the electron liquid. The hydrodynamic approach 
does not assume the Fermi-liquid behavior of the electron 
fluid, but it also applies to classical and semiquantum 
strongly correlated liquids for which experimental real¬ 
ization of the hydrodynamic regime is more realistic. As 
the temperature is lowered, the Fermi-liquid theory be¬ 
comes applicable. At the same time, the equilibration 
length due to intralayer electron-electron collisions be¬ 
comes longer and the system eventually crosses over into 
the collisionless regime of Coulomb drag. 

Previous microscopic calculations of Coulomb drag 
were made under a tacit assumption of the collision¬ 
less regime with respect to intralayer electron-electron 
scattering, namely when the intralayer mean free path 
exceeds interlayer separation £ ^ d. The only excep¬ 
tion is the paper Ref. [45] where crucial importance of 
the intralayer equilibration has been emphasized in the 
context of drag between one-dimensional quantum wires. 
In the two-dimensional case, at large interlayer spac- 
ings, the crossover from the collisionless to the hydro- 
dynamic regime occurs within the range of applicability 
of the Fermi-liquid theory. This enables construction of 
a microscopic description of Coulomb drag in the entire 


crossover interval between the collisionless and the hy¬ 
drodynamic regimes. 

In the present work, we develop a theory of Coulomb 
drag for an arbitrary relation between the interlayer spac¬ 
ing d and the intralayer equilibration length £. We find 
that the crossover from the collisionless, £ ^ d^ to the 
collision-dominated (hydrodynamic) regime, £ d/is hi¬ 
erarchical and is characterized by several parametrically 
distinct energy scales. Depending on the temperature, 
drag is dominated either by the particle-hole continuum 
or plasmon modes. In addition to reproducing the previ¬ 
ously known oc temperature dependence in the colli¬ 
sionless regime [Eq. Q], we clarify the more subtle linear 
in oc T term briefly discussed in Ref. [271 More impor¬ 
tantly, we identify a regime in which the drag resistivity is 
dominated by plasmons and follows a temperature de¬ 
pendence. Apart from being relevant to experiments, our 
work provides an alternative approach to the Coulomb 
drag effect that is based on th e stochastic Boltzmann- 
Langevin kinetic equation.S^^^ 

The paper is organized as follows. In Sec.|n|we present 
a qualitative discussion for the Coulomb drag effect in 
the crossover region from the collisionless to the hydro- 
dynamic regime. In Sec. Ill we derive a general formula 
for drag resistivity using the formalism of the stochastic 
kinetic equation. In Sec. jlVj we focus on the situation 
in which the interlayer momentum relaxation rate due to 
drag, IjTD = is smaller than the equilibration 

rate of the electron liquid, Xjr ^ /Ep^ when the de¬ 
scription of drag is greatly simplified. This situation is re¬ 
alized at sufficiently large intralayer separations. In this 
regime, the electron liquid in the layers may be assumed 
to be in local equilibrium characterized by a drift velocity 
in each layer. To obtain specific predictions for the tem¬ 
perature dependence of the drag resistivity in Sec. |V| we 
adopt a modei coiiision integrai characterized by a singie 
reiaxation rate 1/r. In Sec. VI we anaiyze our resuit in 
various temperature regimes, and we uncover the struc¬ 
ture of the crossover in temperature dependence of the 
drag resistivity. We find that at iow temperatures, drag 
is dominated by the particie-hoie continuum, whereas in 
the coiiision-dominated regime of higher temperatures, 
drag is governed by the piasmon modes. The ratio of 
Pd{T)IT‘^ is a nonmonotonic function of temperature 
that exhibits a broad peak. In Sec. m we discuss our 
findings and compare them with avaiiabie experimentai 
resuits. In our summary, we aiso discuss open questions 
and directions for future research reiated to magnetodrag 
phenomena. 


II. QUALITATIVE DISCUSSION 

Couiomb drag originates from interaction between fluc¬ 
tuations of eiectron densities in the two iayers with spa- 
tiai scaies of the order of the interiayer distance d. At 
iarge interiayer separations, kpd ^ 1, such fluctuations 
invoive a iarge number of particies. To ieading order 
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in the instantaneous electron density at spatial 

scales of order d is given by the mean density of the elec¬ 
tron fluid. In this approximation, however, the Coulomb 
drag vanishes. The fluctuations of electron density that 
result in drag ultimately arise from the discrete nature 
of charge carriers, and at small spatial scales they have 
the character of “shot noise”. At spatial scales of order 
d, the characteristic frequencies and character of propa¬ 
gation of density fluctuations that govern drag, depend 
sensitively on the mean free time of the electrons due to 
intralayer electron-electron scattering, r. 

If the characteristic frequencies of density fluctuations 
responsible for Coulomb drag exceed 1/r, the influence 
of intralayer electron-electron scattering on the dynamics 
of such fluctuations may be neglected. We refer to this 
regime as collisionless. We stress that even in the colli- 
sionless drag regime, electron-electron scattering is essen¬ 
tial for establishing the steady state. Only its influence 
on the propagation of density fluctuations is negligible. 
In the opposite case, i.e., the collision-dominated drag 
regime, electron-electron collisions play a crucial role not 
only in establishing the character of steady-state flow, 
but also in the propagation of density fluctuations re¬ 
sponsible for drag. 

In the collisionless regime, the density fluctuations con¬ 
sist of ballistically propagating particle-hole excitations 
and of plasmons. The characteristic energy scales asso¬ 
ciated with these two types of density fluctuations are, 
respectively. 


Td = i^pi = TdV^. ( 2 ) 

At T <T(i, drag is dominated by the particle-hole contin¬ 
uum and follows the quadratic temperature dependence, 
Eq. Q. The scaling can be simply understood from 
the phase-space argument for the low-temperature two- 
particle scattering arising from the ^ T smearing of both 
the initial and final states near the Fermi surface. The 
power law oc falloff of pD follows from the screening 
properties of the Coulomb potential. 

At higher temperatures, T > drag continues to 
be dominated by the particle-hole continuum but the 
quadratic temperature dependence is replaced by the lin¬ 
ear one, pD because of phase-space limitations 

associated with the predominance of small-angle scatter¬ 
ing. The crossover to the collision-dominated regime for 
the particle-hole excitations happens at 1/r ^ T^. For 
Fermi-liquids this occurs at Tc ^ yjEpTd ^ ujpi. The 
crossover to the collision-dominated regime for plasmons, 
1/r > ujpi^ occurs at a much higher temperature scale 
Th ^ yjEpujpi. Above this temperature the electron liq¬ 
uid enters the hydrodynamic regime where the distinc¬ 
tion between the particle-hole continuum and plasmons 
is no longer meaningful. Drag resistivity in the hydrody¬ 
namic regime followthe 1/T behavior. The crossover 
regime 


ujpi<T < Th 


( 3 ) 


has not been previously explored. Note that the extrap¬ 
olation of the hydrodynamic resulfP^ 


pd{t) ^ ^ / 1 y 

pQ T \kFd) 


( 4 ) 


to temperatures T yields a greater drag resistivity 
than the extrapolation of the ballistic result to by a 
parametrically large factor kpd ^ 1. One concludes then 
that fast intralayer equilibration mediated by electron- 
electron collisions ultimately renders stronger drag effect. 
The detailed comparison of various regimes is given in 
Sec.lVD 

To get some physical insight into the origin of the ap¬ 
parent mismatch between collisionless and hydrodynamic 
answers for the drag resistivity, one should recall that 
hydrodynamic description follows from a more general 
kinetic theory by retaining only zero-modes of the colli¬ 
sion integral. These modes correspond to the conserva¬ 
tion laws of particle number, momentum, and energy of 
the liquid, and they are described by the hydrodynamic 
equations. The more general kinetic description captures 
not only the hydrodynamic modes but also modes that 
are relaxed by electron-electron scattering. Their fluc¬ 
tuations decay rapidly with the typical rate of inelastic 
collisions and govern the crossover regime. 

Below we develop a general theory of Coulomb drag 
in clean double layers with large interlayer spacing that 
is valid in both collisionless and collision-dominated 
regimes, and it enables us to explore the crossover be¬ 
tween them. The dynamics of density fluctuations can 
be described using the g eneralization of the Boltzmann- 
Langevin formalisnPSES to nonideal gases.In this ap¬ 
proach, one describes the state of the system by a fluc¬ 
tuating distribution function f = f ^ 6f averaged over 
a physically microscopic spatial scales (of order d in our 
case) containing a large number of particles. Time evo¬ 
lution is described by the system of equations for the 
average distribution function / and the fluctuations 5f. 
The evolution equation for the fluctuating part Sf is the 
Boltzmann equation with a fluctuating Langevin source 
whose variance is determined by the average distribution 
function /. The short-range part of the Coulomb inter¬ 
action between electrons is described by the collision in¬ 
tegral whereas its long-range part enters the Boltzmann 
equation in the form of an external potential related to 
the density of fluctuation by the Poisson equation. The 
evolution equation for the average distribution function / 
also differs from the standard Boltzmann equation. The 
difference arises because in contrast to the assumption 
of complete relaxation of correlations that underlies the 
Boltzmann equation, only relaxation of short-range cor¬ 
relations is assumed, whereas long-range correlations are 
taken into account. SSI As a result, the evolution equation 
for the average distribution function / contains the corre¬ 
lator of long-range density fluctuations of the distribution 
functions. 
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III. BOLTZMANN-LANGEVIN APPROACH TO 
COULOMB DRAG 

There are three different computational approaches to 
the problem of Coulomb drag (C D) re sistivity, which are 
based on either the Kubo formula,!^^ ^ and closely related 
to it the memory-fu nction fo rmalism,^^^^or alternatively 
the kinetic equationP^HzHH] latter approach has cer¬ 
tain advantages over the other methods since it allows 
one to account on equal footing for inter- and intralayer 
interactions, and in principle it can be generalized to 
nonequilibrium situations. It should be noted that pre¬ 
vious tr eatment s of the CD problem based on the kinetic 
equatiorP^I^^I^ had been carried out in the main kinetic 
approximation, namely without inclusion of the stochas¬ 
tic Langevin forces. Furthermore, the kinetic equation 
for the electron distribution function due to intralayer 
collisions has not been solved explicitly. Instead, the dis¬ 
tribution function was extracted from the argument of 
the Galilean boost in the moving frame with electron liq¬ 
uid and the electron polarization function was then com¬ 
puted for the bare noninteracting limit. Because of these 
reasons, such a description does not adequately capture 
the crossover to the collision-dominated regime. We will 
show that not only is the nonequilibrium part of the elec¬ 
tron distribution function important in the calculation of 
the drag response, but there are also corresponding cor¬ 
rections to the electron polarization function that give 
rise to contributions to drag resistivity. Moreover, the 
functional form of the polarization function is strongly 
affected by intralayer collisions, which thus require sepa¬ 
rate serious considerations. This physical picture has cer¬ 
tain parallels with the problem of nonequilibrium ffuctu- 
ations and the shot noise. Indeed, it is known that elec¬ 
tron collisions strongly modify the spectral function of 
current ffuctuations and change the Fano factor.!^ Drag 
may be interpreted as the rectification effect of nonequi- 
librium ffuctuations of electron liquicP^ so that one nat¬ 
urally expects that equilibration processes play a promi¬ 
nent role for this phenomenon. This observation suggests 
using Boltzmann-Langevin kinetic theory, which was ex¬ 
tremely fruitful in the context of the shot-noise problem, 
in application to the drag effect. 


A. Boltzmann-Langevin kinetic equation 

We consider a Coulomb drag double-layer setup in 
which the active drive-layer is driven out of equilibrium 
by an applied bias while the induced response is mea¬ 
sured in the other passive drag layer. We describe the 
state of the electron liquid within each of the layers la¬ 
beled by an index i by the respective distribution 

functions /* = /* + where /* is the time average 

distribution function and 6f^ is its fiuctuating part. The 
general Boltzmann-Langevin (BL) kinetic equation iJ^ 

[Qt + V . V - e(E^ - V0^)Sp]r = St{r} + 5J\ (5) 


The right-hand side of this equation describes the fiux 
of particles into a given phase-space point due to col¬ 
lisions. The average fiux is described by the colli¬ 
sion integral St{/*}. The stochastic nature of collisions 
causes ffuctuations in the probability fiux, which are de¬ 
scribed by the extraneous source 6J\ The second-order 
correlation function of the extraneous sources is local, 
(JJ(p, r, t)6J{p', r', t')) (X 6{t — — r'), refiecting the 

locality of collisions, and it has been evaluated by Kogan 
and ShulmaiP^ for various types of scattering processes. 
Its specific form depends on the form of the collision in¬ 
tegral, and it can be found in Refs. [47| and |49| 

The fiuctuating electric potentials in Eq. (§ are 
related to the distribution functions in both layers by 
the Poisson equation 

O/jTp r 

= -vr / (6) 

J 

where dV = 2d^p/(27r)^ denotes the density of states in 
two-dimensional momentum space with the factor of 2 ac¬ 
counting for the spin. In addition, in Eq. ^ we also used 
the Eourier representation, and we denoted the interlayer 
distance by d and q = |q|. Isolating the fiuctuating part 
of Eq. (§, we obtain in the Eourier representation 

{-ioj+ivn-eE'‘dp)Sflg+ieq(plgdpf = SSt{f}+SJ^^g, 

( 7 ) 

where (5St{/*} is the fiuctuating part of the collision in¬ 
tegral. To leading order in the ffuctuations, the latter 
may be obtained by expanding the collision integral to 
the linear order in Sp around /b Eor the average part 
of the distribution function, we get from Eq. ^ 

- eE^apf + e{Vc^^d^Sf) = St{f}, (8) 

where angular brackets (• • •) denote averaging over the 
ffuctuations. 

Equations (§, 0 , and 0 supplemented by the ex¬ 
pression for the variance of the Langevin fiuxes SJ'^ 
in terms of the average distribution function /* con¬ 
stitute a closed system describing the ffuctuations and 
Coulomb drag in the double-layer setup. Drag resistivity 
PD = Et/j^ is defined as the ratio between the electric 
field generated in the drag layer in response to the current 
fiow in the drive layer. 

At this point, we briefiy outline the program of calcu¬ 
lations that will lead us to the general expression for pjj. 
Eirst, we need to find the Green’s function of the Boltz¬ 
mann equation Q. Next we will use the Boltzmann- 
Langevin scheme to evaluate the density-density correla¬ 
tions and consequently the dragging force. To perform 
averaging over the stochastic fiuxes, we will assume that 
correlators of Langevin sources have equilibrium form in 
each layer, one moving and one at rest, thus employing 
the fiuctuation-dissipation relation in the boosted refer¬ 
ence frame. 
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B. Drag force 

Using the general approach summarized above, the 
drag resistivity may be determined by solving for the 
nonequilibrium electron distribution function that arises 
in response to a driving (staggered in the layer index) 
electric field and evaluating the resulting staggered cur¬ 
rent. Alternatively, it may be advantageous to evaluate 
the interlayer momentum transfer rate (drag force) for 
a given nonequilibrium distribution of electrons. In the 
present section, we derive a general expression for the 
drag force in terms of the electron distribution function 
in the layers. In Sec. |IV| we use this expression to evaluate 
drag resistivity at large interlayer separations. 

A formal solution of Eq. 0 may be written in the form 

+ < 54 ,,), (9) 


where 

= {-iu: + zvq - eE*5p)(54,, - 5St{f} 

( 10 ) 

is the resolvent of the Boltzmann equation. We note 
that if the stationary distribution function p does not 
correspond to equilibrium, the linearized collision inte¬ 
gral does not generally have zero modes. However, if 
deviations from equilibrium are small, some of the eigen¬ 
values of the linearized collision integral are expected to 
be anomalously small; the latter correspond to hydro- 
dynamic modes. We also note that the source on the 
right-hand side of the above equation is localized in phase 
space in the vicinity of the Fermi surface. 

Multiplying Eq. 0 by p and integrating by parts, we 
obtain 


= e(V 0 Mn*), ( 11 ) 

where we used conservation of momentum in collisions 
f dEpStj/}, and we introduced the average density n* = 
f dTp in layer -2 and its fluctuations Sn^ = f dTSp. 

Using the Poisson equation 0 we can reduce Eq. ( [IT] ) 
to 

= —27re^ ^ - (^6n (g) ^ • (12) 

q,uj 


where we denote the Pauli matrices in the layer space by 
(j. From Eq. 0 we get 

(15) 

Here denotes the extrinsic density fluctuations 

that would be induced by Langevin sources in the Fermi 
gas in the absence of long-range Coulomb interactions 


= J dTpdT,K^^{p,k)Sj:,Jk). (16) 

The polarization operator function 

K,, = J dTpdTkKJp,k)icidpr (17) 


describes the response to a fluctuating potential in a 
nonequilibrium steady state P{p) and also in the pres¬ 
ence of the electric field Eh We can rewrite Eq. (15) in 
the form 


(i - flu;,q%)5n^,q = 


(18) 


where matrices (in the layer index) of polarization oper¬ 
ator and interaction potential are denoted as follows: 

^ ^ 27re^ / 1 p~Qd \ 

n = — (^ gU 1 j • (19) 

From the last two equations, we read off the power spec¬ 
trum of density fluctuations in the form 


{6n0 6ff)uj,q = 

(1 - U^,qVq)-^N^,q{l - (20) 

where Nuj,q = (g) Substituting this ex¬ 

pression into the trace in Eq. (14), we get 


Tr 


= Tr 




-1 


ay 

(y-i - - fi-^,-q)-^ 

By employing Eq. (p^, we notice that 
- 1 ; 






y^q 


\2'Ke^ J 1 — e 


( 21 ) 


( 22 ) 


Here we wrote density fluctuations in the two layers in 
the form of a column vector and introduced the spectra 
power of density fluctuations in the standard way 

{dn^,q®dnl',q') = (27r)^(5(w+w')(5(q+q') {5n®5fF)^^. 

(13) 

Summation over the frequency and momenta implies 
'^q,uj = f ^fd^q- Changing the variables (w,q) -> 
(—cj, —q) under the integral we can rewrite the drag force 
as 


eE*n* = —Tre 


2 qe 


— qd 


q,uj 


-Tr 


eg 


[6n (g) 


(14) 


and therefore the formula for the dragging force reduces 
to 


eE%* 


Tre^ ^ g^eq 
2ehi‘^ ^ sinh{qd) 


Tr 


. /g._i 






where eq is the unit vector in the direction of q. Here 
we also introduced the dimensionless interaction matrix 
'dq = (6g/27re^)Ug, the dimensionless polarization opera¬ 
tor P = z/“^II, where u is the single-particle density of 
states, and the inverse Thomas-Fermi screening radius 
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n = 2'Kve^/e. The above equation expresses the drag 
force in terms of the density response functions in each 
layer, and the correlator of extrinsic density fluctuations. 
We note here that the latter in general can not be re¬ 
duced to a density response function. In equilibrium, 
and are proportional to the unity matrix, and 
the interaction matrix Vq has only do and components. 
As a consequence, the trace in the integrand of Eq. (23) 
vanishes resulting in zero drag force, as it should be. 

For the further convenience with the intermediate steps 
of calculation we introduce the following matrices 




qe' 


qd 


2k, sinh(g(i) 




■ <^5^’ 


(24a) 


= 


qe 


qd 


2k sinh(gd) 


^0 P—u!, — q’> 


(24b) 


V = 


Q 

2Ksmh.{qd) 


(24c) 


We remind the reader that the polarization operator ma¬ 
trix is diagonal in the layer index, however it is propor¬ 
tional to do only for the case of identical layers. We will 
assume this case for the simplicity of further considera¬ 
tions, and we will provide the necessary generalization for 
the case of unequal layers at the end of this section. With 

these new notations we observe that — P^^qj = 


and (^qi - P_„,_g) ' = {P*-V)-\ Next, 

let us denote nonequilibrium corrections to various quan¬ 
tities by A so that to the linear order in those we get from 

(d) 


eE%* = 


ire 




2eK^ ^ sinh(gd) 


(p* _ V)-^cfy{P - V)- 


AP(P - - V)“^AP*}] . 


(25) 


may be characterized by two hydrodynamic velocities in 
each layer, and deviation from true thermal equilibrium 
is characterized by the difference between the layer ve¬ 
locities. 


IV. DRAG AT FAST INTRALAYER 
EQUILIBRATION 


Drag resistivity may be characterized by the relax ation 
rate Ifro = of the staggered momentum.l^^Ml 

At large interlayer separations, this rate is significantly 
smaller than the intralayer equilibration rate of the elec¬ 
tron fluid, 1/r, ^ pDTie^/m. Under this condition, 

the electron distribution function in each layer is well 
approximated by the equilibrium distribution with the 
drift velocity ±u/2 corresponding to the current density 
in the layer. In this case, we can obtain the nonequi¬ 
librium parts of various quantities by applying Galilean 
boost, namely 

A-FL,? = ^^zduiPoj,q, AP<^,q = (26) 


where Puj,q is the matrix distribution function of collec¬ 
tive bosonic excitations, which at equilibrium is given by 
= coth(cj/2T). Since by the fluctuation-dissipation 
theorem 


^uJ,q — 2^ 


P — P 

uj,q —uj, — q 


B,. 


(27) 


we readily find in the notations of Eq. ( |24a[ ) 

{aP-aP*)b + {P-P*)abY (28) 


= I 


Using this result in Eq. (25), and after some matrix al¬ 
gebra, we can rewrite the trace as a sum of two contri¬ 
butions. 


eE^n* = 


Ea77^['^i+'^ 2], (29a) 


I'Kve 
4e/^^ smh{qd) 


This equation constitutes the essential result of this sec¬ 
tion. It relates the driving electric field to the polariza¬ 
tion operators of the individual layers and the correlators 
of Langevin fluxes in them. The latter are determined by 
nonequilibrium distribution functions, which are in turn 
determined by the driving electric field via the system of 
equations (§-(§. Further calculations of drag require 
explicit forms of the nonequilibrium corrections to the 
polarization function and the correlator of the Langevin 
fluxes. Determination of the nonequilibrium distribution 
function in the general situation is a very difficult prob¬ 
lem. The situation simplifies dramatically when the rate 
of momentum transfer between the layers due to drag is 
smaller than the equilibration rate of the electron liquid 
due to intralayer electron-electron scattering. This is al¬ 
ways the case when the layers are sufficiently far apart. 
In this situation, the nonequilibrium state of the system 


Ti = B^Tt 


ay[A{P-v)-^ 


A(p*-V)-i] , 


(29b) 


- vyH'P - - V)“^ . 

(29c) 

It is easy to see that the variations A{P — V)“^ and 
A(P* — V)“^ contain only matrices do, dx and dz- There¬ 
fore, the trace of their product with dy must vanish iden¬ 
tically, Ti = 0. The second trace does not vanish because 
of the explicit presence of the dz matrix arising from the 
staggered boost in the layers. Indeed, after cyclic permu¬ 
tation of matrices under the trace, we observe that 




(P* _ V)-i,Tj,(P - V)-1 = 

ay{\r\^-V^) + iaMP-P*) 


|P-V|2|P +V|2 


( 30 ) 
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where in accordance with Eqs. (24a) and (24c) we wrote 
V = 7^(Jo and V = with V = — P. As a conse¬ 

quence of the above equality, we find for the nonvanishing 
contribution to the trace in the formula for the drag force 


To = 


z(qu) 


V{T-T*f 


2Tsmh^{uj/2T) \V - V\^\V + Vp 


(31) 


At this point we are prepared to find drag resistivity. 
Since averaging over the fluctuations has been performed, 
all the quantities appearing in the expression for the drag 
force should be understood as being equilibrium aver¬ 
ages. In particular, to simplify notations in what follows 
we relabel the time average density n ^ n with the un¬ 
derstanding that n is the equilibrium density of carriers 
in quantum wells. 


We collect all the factors in Eqs. (29) and (31), notice 


that upon angular averaging over the momentum trans¬ 
fer (eq(equ))q = u/2, aud thus we obtain for the drag 
resistivity po = E/enu the following result 


PD 

PQ 


1 


E 




{ImVf 


S'kTv? sinh^(w/2T) \'P — Vp|P ■ 


■vp 


(32) 


This expression can be brought to a more familiar form 
by observing that IxaP = —ImP, and consequently 


V(p_p*)2 
|P-V|2|P +V|2 “ g3 

= _^e-«‘^Im ( A) Im A 

4 q \S+J \S. 


(1 - e-2«‘i)e-«‘^(P - P*)2 




(33) 


where in a usual way we introduced the dielectric func¬ 
tions corresponding to the symmetric and antisymmetric 
density modes 

£:±(a;,g) = l-(34) 


As a result, Eq. (|32|) can be equivalently rewritten as 
1 


PD ^ _ 

PQ ~ 2567rTn2 


^ q^e ‘^^Im(£:^^)Im(£:_^) 


q,uj 


sinm(u;/2T) sinh(g'(i) 


It is instructive to make several comments in regards 
to Eq. ( [ 3 ^ . {%) If one neglects the plasmon resonances 
in the denominator of Eq. (32), then one finds that drag 
is simply proportional to (ImP)^, which is a standard 
formula that applies to the contribution of the particle- 
hole continuum, (n) In general, the full expression (33) 
gives a complete formula that accounts for plasma res¬ 
onance of the two dielectric functions. As will be ex¬ 
plained later, plasma modes govern drag resistivity in the 
collision-dominated transport regime, {in) In the case of 
non-identical layers, Eq. (32) should be generalized as 
follows: in the numerator, one should replace Im(P)^ 
ImP^ImP'*^ and ^ while in the denominator, 

v? n^ri^ and |P —+ ^ |Det(P —V)p, where 

the determinant can no longer be simply factorized as a 


product of two independent terms for each of the layers. 
In our analysis we concentrate on the case of symmetric 
layers. Even though the analytical calculations are still 
doable for the general nonsymmetric case, the obtained 
results for the drag resistivity become quite cumbersome 
and do not bring any new physics insight concerning the 
temperature dependence of the transresistivity. 


V. RELAXATION TIME APPROXIMATION 


When analyzing temperature dependence of the drag 
resistivity from Eq. (32) we must discuss analytical struc¬ 
ture of the polarization function in (cj, g')-plane. In 
most of the previous studies this function was calcu¬ 
lated for the bare noninteracting limit with respect to in- 
tralayer electron-electron collisions. This approximation 
is only sufficient to describe the low-temperature colli¬ 
sionless regime where drag is dominated by the particle- 
hole continuum. The finite intralayer equilibration length 
^ changes in a significant way, which becomes quan¬ 
titatively important for drag already at moderately high 
temperatures T > vp/d. Technically, finite i makes 
ImP^^g 0 8it the high-frequency limit, uj > vpQ^ where 
the conventional collisionless result the yields vanish¬ 
ing spectral weight of particle-hole excitations. Conse¬ 
quently, as a result of equilibration processes, plasmon 
poles lie within the high-frequency tail of the an particle- 
hole continuum of states and yield important contribu¬ 
tion to drag resistivity. This mechanism becomes pro¬ 
gressively stronger as temperature gets higher, with thus 
a shorter equilibration length, and eventually plasmon 
modes take over and dominate the Coulomb drag. These 
points will be explicitly illustrated by the subsequent cal¬ 
culation, while we turn our attention now to the deriva¬ 
tion of P^^q at finite £. 

We treat the intralayer collision integral in the relax¬ 
ation time approximation 


St{/} = --{Sf - fh)- 


(36) 


We assume that r is dominated by electron-electron in¬ 
teractions, so that (X T^. Importantly, the Eermi- 
liquid-like quadratic temperature dependence of the re¬ 
laxation time has been experimentally confirmed even in 
the regime of correlated electrons with thus strong inter¬ 
actions Tg ^ 1.1^ In Eq. (36) fh stands for the hydrody¬ 


namic part of the nonequilibrium distribution, which has 
the locally equilibrium form 


fh = 


6p + m(u • v) + 


ST{s - p) 


dsfo (37) 


where /o = + 1]“^ is the Eermi distribution. 

The Lagrange multipliers 6p^ u and 6T have the mean¬ 
ing of the local change in, respectively, the chemical po¬ 
tential, hydrodynamic velocity, and temperature. Their 
values are found from the conditions that the collision 






















integral (36) conserves the particle number, momentum, 
and energy: 


J dTSf = -5^1 j dTdefo = vdji, (38a) 


drp(5/ = - 


J dTp{pu)de 


fo = mnu, 


(38b) 


I dris-^)Sf = -^-^ J dr{e-^) 


(38c) 

with u = d/^n being density of states. With the collision 
integral in the form of Eq. (36) the solution for the lin¬ 
earized Boltzmann-Langevin kinetic equation 0 can be 
presented in the form 




q . . I ' ^e/o 

1 — lUJT -1- irqy V 


rSJuj 


ieT(l)q^^qv + Sfi^^q m^ru^^q -f 


{s - IJ)8T^ 


(39) 


where we suppressed the layer index-z for brevity and 
used the Fourier representation. We substitute the above 
expression into the integral constraints of Eqs. (38). It 


is convenient to separate the integration over the phase 
space dr into the integration over the energy and the an¬ 
gular averaging over the Fermi surface, (...). The latter 
can be performed with the aid of the following relations 


1 


1 — iujTiTC\y/ D^ q^ 

V \_ _ _^q 

1 — iujT + irqv 




1 - 


1 — iujT 


QiQj 


(40a) 

(40b) 


1 MLu.,q-M^^q)^^M^Jij,{AOc) 

1 — + 2rqv/ 


where we introduced the following dimensionless quanti¬ 
ties 


D^^q = y/(l - iur)‘^ + ('q^a) 

^ “ (1 “ (41b) 

-^<^.9 = ^^(-^“.9 “ 1 (41c) 

When performing integration in dT over the absolute 
value of V we make use of the advantage that in the de¬ 
generate electron liquid all integrals are dominated by 
momenta near the Fermi surface. As a result, we ob¬ 
tain a closed system of linear algebraic equations for the 


Fourier components of Lagrange multipliers 


I 1 - - 1 O/J^^^q ^ IV [I -' ^ 




uj,q 


n 


uj,q 


q^T 


^ . 1 - iujT , 

^^4^00,q ( 1 ^ I T 


h 


TS JqdX^ 

— iuiT + irqv ’ 


P \ I - - iT{n - vmC^^q)cpi^ q = 


n 


uj,q J 


2 2//. f 

eVT q (Puj,q^uj,q / Z T 

J 1-t 


iT^vqSJ^^qdT 




iujT + irqv 

T{e - iJ.)6J^^qdr 
— IUJT + irqv ’ 


/? 


(42) 

(43) 

(44) 


where V^^^q and C^^q are obtained from D^^q and L^^q 
with the replacement of v by the Fermi velocity v ^ vp. 
It is interesting to observe that fluctuations of tempera¬ 
ture decouple from the fluctuations of the density and 
drift velocity. Since our goal is to get the polariza¬ 
tion function for the density fluctuations, we only need 
the first two of the above equations. Notice that the 
latter are also consistent with the continuity equation 
—iujdn^^q + inqu^ g = 0. By excluding u^^^q in favor 
of dn^^q = vdfji^^^q from the above system, we can cast 
the equation for the density fluctuations in the form of 
Eq. (15) and consequently read off the polarization oper¬ 
ator 



(45) 

This important result together with Eq. ( [3^ enables us 
to study the temperature dependence of the drag resistiv¬ 
ity in various transport regimes, which is our immediate 
goal. 


VI. DRAG RESISTIVITY 

We have determined that there are three distinct tem¬ 
perature scales, and consequently four separate regions, 
where Coulomb drag resistivity exhibits qualitatively dif¬ 
ferent behavior. At lowest temperatures, the equilibra¬ 
tion length is very large, i = vpT ^ d, electron kinetics 
is essentially collisionless, and drag is dominated by the 
particle-hole excitations with characteristic wave vectors 
q ^ d~^ and frequencies uj ^ T. Since the spectral edge 
of the particle-hole excitations is set by the line uj = vpq^ 
this introduces the characteristic crossover temperature 
Td ^ Ep/kpd [see Eq. 0 ]- 

As temperature increases, the equilibration length gets 
shorter and may become comparable to the interlayer 
spacing, I ^ d^ which happens at the characteristic tem¬ 
perature 


Tc ^ Ep^Jkp /ti^d. 


(46) 
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This scale signifies the beginning of the collision- 
dominated transport regime. The remaining scale is de¬ 
duced from the condition when the equilibration rate 
becomes comparable to the energy scale of plasmon 
modes, ujpi ^ at characteristic for the drag wave 

vector q ^ d~^. This yields the temperature scale 

(47) 


above which hydrodynamic transport regime sets in. Our 
goal is to describe the crossover in the temperature de¬ 
pendence of drag resistivity throughout the entire range 
of hierarchical scales, < Tc < Before we pro¬ 
ceed with task it will be convenient to introduce sev¬ 
eral dimensionless variables and parameters. We will be 
measuring momenta q and frequencies uj in units d~^ 
and Td = vp/d respectively, thus introducing x = qd 
and w = uj/Td. We also introduce dimensionless equi¬ 
libration length I = vpT jd^ and dimensionless frequen¬ 
cies for symmetric (antisymmetric) 

plasmon modes. In these notations equation ( [3^ can be 
rewritten in a manifestly dimensionless form. 


PD 

PQ 


Td 1 X 

327rT {nd‘^)‘^(tvd)^ sinh^(x) sinh^(rcT^/2T) 




p 


-1 

W,X 


x‘^ P 


-1 

w,x 


(48) 


We need to explore now various asymptotic limits of this 
formula. 


A. Collisionless regime T < Tc 


In the low temperature limit drag resistivity is domi¬ 
nated by the low-energy excitations since plasmon poles 
are much higher in energy. Indeed, for typical x < 1 we 
have Ws(^a) ^ 1 while integral in Eq. (48) is dominated 
by frequencies w ^ '?^s(a )5 so that we can replace the 
denominator in the second line of Eq. (48) to unity and 
obtain a simplified expression 




D 


1 


E 




Pq 327rT {nd?Y{i<idY sinh^(x) sinh^(rcT^/2T) 

^(49) 

In the current notation with the superscript-ph in we 
want to emphasize that this contribution originates from 
particle-hole modes to distinguish it from contributions 
due to plasmons that will be discussed below. 

Eor T Tc equilibration length is / ^ 1, and polar¬ 
ization operator in Eq. (Eol can be taken in the main 
approximation. Erom Eq. (|45p we find in the dimension¬ 
less variables 


P — 1 — 

J- W.X - 


W 


— x‘^ 


(50) 


Eor this temperature regime, keeping correction terms to 
Pw^x in powers of 1// will result only in subleading cor¬ 
rections to drag resistivity in Eq. (49) in a small param¬ 
eter T/Ep ^ 1. At the lowest temperatures, T < Td^ 
frequency integration in Eq. (49) is dominated by the 
range where w x, so that it is sufficient to take 
Im(P“^) ^ w/x, and thus we obtain 


pt 


7>'C(3) 


PQ 16 {kpciyiKciy EP 


(51) 


which is a well-known result 

At the higher temperatures, T > Td^ we can approxi¬ 
mate sinh^(rcT^/2T) « {wTdl2T)‘^^ however we need now 
the full expression for Pw,x from Eq. (50) since the lead¬ 
ing contribution to integrals comes from w ^ x. As a 
result, drag resistivity crosses over to the linear temper¬ 
ature dependence 




1 


Pq 360 {kpd)^{hzd)‘^ Ep 


(52) 


It is worthwhile emphasizing that the dependence on the 
interlayer separation changes as well from being propor¬ 
tional d~^ to d~^. 

To account for the contribution from plasmon reso¬ 
nances to drag resistivity in this temperature regime, we 
need a more accurate expression for the polarization func¬ 
tion in the high frequency limit. Expanding Eq. (45) 
under the assumption that x/{w Pi/I) <C 1 one finds 


P~^ = 


2w‘^ 

x‘^ 


- I 


w 


2{w + i/l) 


(53) 


This expression allows one to approximate resonant de¬ 


nominators in Eq. (48) by plasmon poles 


- P~^ - 1 

2wl 


-2 




[{w - Waf + ll][{w + Waf + 7a] 

(54) 

where we have introduced plasmon damping rates 


I 


la = 


8 1 + {wjy ’ 


a = s,a. 


(55) 


The expression (55) is based on the relaxation-time ap¬ 


proximation and describes the plasmon attenuation rate 
in the entire crossover between the hydrodynamic {wq^I 
1) and collisionless {wo,l ^ 1) regimes. In the collisionless 
regime, it is consistent with the results of the microscopic 
treatment of plasmon attenuation rates in Ref. [52l 


Eor T <^Tc the plasmon contribution to pp in Eq. (48) 
is dominated by the pole at the frequency of the antisym¬ 
metric plasmon. Using Eqs. ([^ and (54) in Eq. (48) 


pi 

Pd. 


Td 


E 


Pq 1287rT {nd‘^Y{i<idY ^ sinh^(x) sinh^(rcT^/2T) 


(wp 


Wn 


[1 {wiPY [{w - WaY + 7a] [(^ + '^aY + ll. 


■•(56) 
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Integration over the frequency here is straightforward 
since poles near ±Wa are narrow and each Lorentzian 
function is effectively a delta-function with the weight 
7r/7a. The remaining x-integration is dominated by the 
range where x <1 where one can approximate w‘^ = hzdx‘^ 
and find as a result 


pg 3C(3) 1 / T y 

PQ 16 \Ef J 


(57) 


This contribution is obviously smaller than that given by 
Eq. (52) in the same range of temperatures. We con¬ 
sidered this term here in such details not only for com¬ 
pleteness of our analysis. It will be shown in the next 
subsection that Eq. (57) will cross over into a different 
power-law behavior of drag resistivity that will overcome 
the respective contribution of the particle-hole contin¬ 
uum. 


B. Collision-dominated regime T > Tc 


Erequency integration in this formula can be done ex¬ 
actly, however it leads to an extremely cumbersome ex¬ 
pression. A much more useful result can be obtained 
by exploring the following simplifying observation. The 
relevant wave numbers for the integrand function in the 
above expression are logarithmically large, x r^ln{Fd/P). 
In this case, the frequencies of the symmetric and anti¬ 
symmetric plasmons nearly coincide, and consequently 
their dumping rates become almost identical. Thus one 
can justify replacing 7^ ^ 7 = (x^//8)[l + tvdl‘^x/2]~^ 
and also P[1 + (87/x^)^. We integrate 

over w by leading poles and make use of the follow¬ 
ing algebraic identities: {wsWa)‘^ = sinh(x), 

{ws^Wa)‘^ = K,dx{l ± \/l — w‘^-\-Wg = K,dx^ which 

leads us to the following intermediate result: 


E 


x^e^ 


_L 

pD 327TTci {nd?‘Y{i<id) ^ sinh(x) 1 + tidl‘^xl2 


p3g2a:^2(i ^ — e 

64/^(i(l + i<idl‘^x/2y 


(61) 


We proceed to consider the intermediate and high tem¬ 
perature intervals. In this case the particle-hole contin¬ 
uum and the plasmon contributions are described by well 
separated peaks in the integrand of Eq. (32) and can be 
considered separately. 

The particle-hole contribution is dominated by the fre¬ 
quency range where re > x so that we need to use Eq. (53) 
for the polarization operator. Eurthermore, we may also 
approximate sinh(reT^/2T) ^ wTdl2T and set plasmon 
denominators in Eq. ([3^ to unity, which thus yields 


It is useful to observe at this point that FdP ^ (T/^/T)^ 
so that the remaining integral can be estimated for two 
asymptotic regions. At the range of intermediate tem¬ 
peratures, Tc <^T ^ T/j,, we estimate 

^ -J—H urn 

f,Q~(kv<lfEy * 

whereas in the hydrodynamic regime of high tempera¬ 
tures T ^ Th we estimate with logarithmic accuracy 




1 


E: 




PQ 2,2'KTd {nd?Y{KdY ^ sinh2(a;) [1 + (^0^]^' 

(58) 

The remaining momentum and frequency integrations 
are elementary here, and one finds 


f^D _ 3C(3) 


1 


Ep 


PQ 16 {kpdfindY T 


(59) 


It is easy to check that at T ^ Tc this result matches 
with Eq. (52), and at the same time it shows that the 
contribution of the particle-hole excitations to the drag 
resistivity becomes a decreasing function of temperature 
in the collision-dominated regime. 

The situation with plasmon contributions is different. 
Accounting for both symmetric and antisymmetric plas¬ 
mon resonances, and using Eqs. (53) and (54), we obtain 
from Eq. (32) 


/og _ T 1 ^ a;2 P 

^ “ 327rTd {nd^Y{KdY ^ sinh2(a;) [1 + 

X TT __ f60) 


pi 

Pd 


Ep 


PQ (kpdPindP T 


(63) 


in agreement with Ref. |44] [see the comment in Ref. [53] . 
We can check that the low-temperature plasmon con¬ 
tribution ( [57| ) matches with Eq. (62) at the expected 
crossover scale, T ^ Tc, while Eqs. (62) and (63) are 
of the same order at T ^ Th where drag resistivity 
reaches its absolute maximum. It is instructive to com¬ 


pare particle-hole [Eq. ( |5^ ] and plasmon [Eq. (62)] con¬ 
tributions to drag resistivity in the collision-dominated 
regime. In particular, at T ^ one easily finds 
that plasmons dominate by a parametrically large fac¬ 
tor p^j^/p^j^ Fd'> 1. 


VII. DISCUSSION 

We have developed a general computational scheme 
to describe nonlocal transport in interactively cou¬ 
pled double-layer systems. As an alternative route to 
existing formulations, our approach is based on the 
stochastic Boltzmann-Langevin kinetic theory. Using 
this approach, we reproduced the known results for 
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Pd 

Pmax plasmons 


P2 


particle-hole continuum 

Pi T 

P 



1/T 


Tc T, 


FIG. 1. Schematic illustration of the temperature depen¬ 
dence of drag resistivity in various regimes. The four intervals 
from low to high temperatures are described by Eqs. (51), 
( [ 5 ^ , ( [6^ , and ( [6^ , respectively. At the crossover scales 
Td and Tc^ the estimated magnitude of drag resistivity is 
Pi/PQ ~ l/{kFd)^{K,d)‘^ and P2IPQ ~ kp / d / (kpdY{ndY , 
while the maximum value of drag, pmax is estimated in 
Eq. (1^. 


the drag resistivity of clean double layers in the low- 
temperature regime, and we studied the previously unex¬ 
plored collision-dominated transport regimes. Our main 
analytical result is Eq. (35) for the drag resistivity, which 


is expressed in terms of the momentum- and frequency- 
dependent dielectric function of individual layers. This 
expression has some parallels with the Lifshitz theory of 
van der Waals forces.^ It holds as long as the intralayer 
equilibration rate is faster than the interlayer momen¬ 
tum relaxation rate. Note that this assumption is valid 
at large interlayer distances, kpd ^ 1. Under these con¬ 
ditions, the electron distribution in the driven system 
is given by the locally equilibrium distribution function 
in each layer and is characterized by the corresponding 
drift velocity. The correlators of Langevin forces driv¬ 
ing the density fluctuations can then be extracted from 
the fluctuation-dissipation theorem. Under more com¬ 
plicated nonequilibrium conditions, drag should be com¬ 
puted from the more general Eq. ( [2^ while the distri¬ 
bution function and Langevin forces should be found by 
solving the corresponding kinetic problem. It is worth 
stressing that the characteristic time scale for electron- 
electron scattering Tee can be significantly shorter than 


that for the impurity. 


or phonon, scattering. 


For example, for relatively clean GaAs quantum wells 
at the doping level ^ and at temperatures 

^ lOK electron-electron scattering time is estimated to 
be Tee ^ 10“^^s whereas and Tph are typically of order 
the order of ^ 10“^s. The relation Tee ^ jus¬ 

tifies applicability of our result for drag resistivity (l3^. 


The general expression (35) applies not only to Fermi 


liquids but also to classical and semi-quantum strongly 
correlated liquids. To make further progress we focused 
on the Fermi-liquid regime. To obtain concrete expres¬ 


sions for the drag resistivity we treated electron-electron 
collisions in the relaxation time approximation. As a re¬ 
sult we expressed the drag resistivity in terms of the sin¬ 
gle relaxation rate, 1/Tee ^ /Ep. We found that there 
are four distinct temperature regimes in which drag ex¬ 
hibits a qualitatively different behavior and is governed 
by different types of density fluctuations - quasiparticles 
and plasmons. The main findings of our work are sum¬ 
marized in Fig. ([^, which schematically illustrates the 
temperature dependences of drag resistivity in the vari¬ 
ous transport regimes. 

At low temperatures, T <Te [see Eq. (§]> drag is dom¬ 
inated by the continuum of quasiparticle excitations with 
collisionless dynamics. Within that temperature domain, 
pD crosses over from being quadratic in T below Td to 
a linear dependence in T above T^. Analytical expres¬ 
sions for these regimes are given by Eqs. (51) and (52), 
respectively. Interestingly, it has been known from the 
early measurements that Pd/T^ ceases to be a mono¬ 
tonic function of temperature already at low tempera¬ 
tures. Indeed, it has been pointed out in Ref. [4] that 
the ratio Pd/T^ displays a noticeable falloff above a cer¬ 
tain crossover temperature. Furthermore, it has been 
pointed out that the ratio d^pp/T"^ appreciably depends 
on d above that temperature, and that the crossover scale 
itself shifts to lower temperatures for the larger inter¬ 
layer separations. All these features are accounted for by 
our results for particle-hole contributions to drag with 
the identification that expected crossover takes place at 
T ^ Ep/{kpd?). 

The collision-dominated regime for quasiparticle dy¬ 
namics sets in at T > Tc- In this regime. Coulomb drag 
is dominated by the plasmon contribution. The temper¬ 
ature dependence of this contribution is very sensitive 
to the plasmon attenuation rate. In the relaxation-time 
approximation, the latter is given by Eq. (55), which is 


consistent with the result of a microscopic calculation of 
Ref. |52j It is worth noting that in a wide temperature 
interval, Te < T < [see Eqs. (46) and (47)], the fre¬ 
quency of the plasmon modes contributing to drag, Upi 
in Eq. (§ , still exceeds the rate of electron collisions, so 
that their hydrodynamic treatment is inapplicable. In 
this intermediate regime, drag resistivity is oc T^, [see 
Eq. (62)]. At T > T/i, the rate of equilibration exceeds 
Upp In this regime, drag can be treated using the hy¬ 
drodynamic approacfP^ and follows the 1/T temperature 
dependence [see Eq. (63)]. In the entire range, drag may 
be represented by the interpolation formula 


Pd 

PQ 


1 


(64) 


cdkpdYpdfT/EF + C2{kFdYE%/T ’ 

with numerical coefficients of the order of unity, Ci ^ 
C 2 ^ 1. This expression allows us to estimate the maxi¬ 
mum value of drag resistivity. 


pmax 

Pd 


Y K/kF'sfKd 


(65) 


PQ pdYikpdY' 

We note that the measurements of Refs. [221 and SOI 
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in the low-density strongly correlated samples revealed 
that within the broad experimentally probed tempera¬ 
ture regime, Pd oc where power exponents vary 

in the range 1 < ar < 3 and 2 < < 5. Perhaps even 

more importantly, the observed magnitude of drag resis¬ 
tivity was one to two orders of magnitude larger than 
expected on the basis of a simple extrapolation of the 
small Vs results. Such a manifestly nonquadratic tem¬ 
perature dependence of drag resistance combined with 
the unexpected magnitude of the effect was attributed 
in the literature to a possible non-Fermi-liquid behavior 
of strongly correlated liquids. In contrast, our conclu¬ 
sions are different. We have found that fast equilibration 
promotes a stronger drag effect, and at the same time 
it is responsible for the nonmonotonic temperature de¬ 
pendence of drag, so that most of the observations find 
a natural explanation within our theory. A very pecu¬ 
liar detail, which should be stressed once again, is that 
counterintuitively plasmons start to dominate the drag 


even at moderately high temperatures, and definitely at 
temperatures way below the plasmon resonances. 
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